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Abstract. Quark mass corrections to the r hadronic width play a significant role only for the strange 
quark, hence providing a method for determining its mass. The experimental input is the vector plus 
axial- vector strange spectral function derived from a complete study of tau decays into strange hadronic 
final states performed by ALEPH. New results on strange decay modes from other experiments are also 
incorporated. The present analysis determines the strange quark mass at the Mr mass scale using moments 
of the spectral function. Justified theoretical constraints are applied to the nonperturbative components 
and careful attention is paid to the treatment of the perturbative expansions of the moments which exhibit 
convergence problems. The resuh obtained, ms(M^) = (120 ± llcxp ± 8v„, ± 19th) MeV = (120l26) ^^V, 
is stable over the scale from M-r down to about 1.4 GeV. Evolving this result to customary scales yields 



m,(l GeV^) = (160+;°) MeV and m,(4 GeV") = (116 



+20-, 
25) 



MeV. 



1 Introduction 

The high precision data on r decays collected at LEP 
and CESR provide a powerful tool to investigate the dy- 
namics of strong interaction at low energies and to deter- 
mine basic parameters of the theory. The QCD analysis of 
the nonstrange inclusive t decay width has 
led to accurate measurements |Q,^|l^,|l2| of the strong 
coupling constant at the r mass scale, as(M^), which com- 
plements and competes in accuracy with the high preci- 
sion determination of as{M'^) from measurements of the 
Z width at LEP. 

More recently, experimental studies of the Cabibbo- 
suppressed width of the r became available |l3[0, al- 
lowing to initiate a systematic study of the corrections 
induced by the strange quark mass in the r decay width 
|,^|l|,|l3,|l|,[l9[|2^,|l|,||. From the separate measure- 
ment of the strangeness S = and S = —1 t decay 
widths^ it is possible to pin down the SU(3) breaking ef- 
fects and to perform a reliable determination of the strange 
quark mass. 



The ALEPH data 1 13 have been used in previous anal- 
yses to extract the value of ms(M^), leading in some cases 
to different, albeit not inconsistent results |]l3|, p^|20| , 22 1. 
In this paper we present an updated common analysis |23|, 

" Now at TRIUMF, 4OO4 Wesbrook Mall, Vancouver, B.C., 
Canada V6T 2 A3 

^ Throughout this paper, charge conjugate states are im- 
plied. 



ta king into account recent experimental information [ p4[ 
^,^6[ in addition to the ALEPH data. Particular atten- 
tion is given to the analysis of theoretical uncertainties 
and to the stability of the results. 

Even within the relatively large statistical errors of the 
present data, the strange quark mass determination from 
T decays has already achieved an accuracy good enough 
to substantially reduce the range quoted by the Particle 



Data Group |24 



2 Theoretical framework 

The theoretical analysis of the hadronic r decay width 
involves the two-point correlation functions 



i7^.>(g) ^^jd'x e'^- (0|T {Vt,{x) V.^iO)^] |0) , (1) 
77^.;'^ (g) ^ ^ / A'^xe^^^ (0|T {A^l^{x) Aj;.(0)t} |0) , (2) 



associated with the time-ordered vector, Vl^^{x) = qj'j'^qi, 
and axial- vector, A^j{x) = qjj'^jzqi, colour-singlet quark 
currents. The subscripts i, j denote the corresponding light 
quark flavours (up, down, and strange). The correlators 
admit the Lorentz decompositions 

+ q>^q'^ ny,j,{q^), (3) 
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where the superscript in the transverse and longitudinal 
components denotes the corresponding spin, ,7=1 (T) 
and J — (L), in the hadronic rest frame. 
The hadronic decay rate of the r lepton, 



Rr 



Ft- 



hadrons t^r(7)) 



Rt.V + Rt,A + Rt,S , 



(4) 



can be expressed as an integral of the spectral functions 
Im7T^(s) and Im7T^(s) over the invariant mass s of the 
final-state hadrons, where 



J ( „\ — 



VudV 

f \Vus 



n. 



(5) 



with \ Vij \ the corresponding Cabibbo-Kobayashi-Maskawa 
(CKM) quark mixing factors. In Eq. (^), Rt.v and Rt A 
correspond to the two terms proportional to in (||) 

and Rt.s contains the remaining Cabibbo-suppressed con- 
tributions. 

Using the analytic properties of U'-' (s), one can express 
Rt as a contour integral in the complex s-plane running 
counter-clockwise around the circle \s\ 



Ml 



Rt 



-TTl 



1 



ds 



L+T 



is)+iD-is)'> . (6) 



We have used integration by parts to rewrite Rt in terms 
of the logarithmic derivatives of the relevant correlators. 



(7) 

(8) 



which satisfy homogeneous renormalization group equa- 
tions (RGE). This eliminates unwanted renormalization- 
scheme dependent subtraction constants which do not con- 
tribute to any physical observable. 

For large enough — s, the contributions to D'^ (s) can be 
organized using the Operator Product Expansion (OPE) 
in a series of local gauge-invariant scalar operators of in- 
creasing dimension D = 2n, times the appropriate inverse 
powers of — s. This expansion is expected to be well be- 
haved along the complex contour \s\ = , except for the 
crossing point with the positive real axis |M . As shown in 
Eq. (|^), the region near the physical cut is strongly sup- 
pressed by a zero of order three at s = M^. Therefore, the 
uncertainties associated with the use of the OPE near the 
time-like axis are expected to be small. Inserting the OPE 
in Eq. (^) and evaluating the contour integral, one can 
rewrite Rt as an expansion in inverse powers of Q. 
This leads to a rigorous prediction for Rt and its different 
flavour components, which, in particular, allows as(M^) 
to be accurately determined [pi[lC|,[ll|,[l2|. 



The measurement of the invariant mass distribution of 
the final state hadrons provides additional information on 
the QCD dynamics. The moments m 



Rt^ = 



Mt 



ds \l- 



M^J \M^J ds 



(9) 



can be calculated theoretically in the same way as Rt = 
The corresponding contour integrals take the form 

J\x\ = l ^ 

+ iTt\x)D^{M^x)}, (10) 

where all kinematical factors have been absorbed into the 
kernels J-^''_^_rp{x) and !F^{x). Their explicit expressions 
were given in Ref. |pj| . Performing the contour integra- 
tion, the result can be written as 



R 



kl _ 



3 [iKdP + IK.P] ^Ew 1 + <5ew + '5'°' (11) 



E ( 



cos 6c 6, 



kl{D) 



sin^^c 



.kl[D] 



r>=2.4, 



where sin^^c = |14sP/[|KidP + iKisP] and where we have 
pulled out the electroweak radiative corrections Sew = 
1.0194 ±0.0040 H and ^J^^ - 0-0010 H- 

The dimension-zero contribution 5^"^ is the purely per- 
turbative correction neglecting quark masses, which, ow- 
ing to chiral symmetry, is identical for the vector and 
axial-vector parts. The symbols ^^j'*-^' = [^^j' y*^ +<5fj"^^]/2 
stand for the average of the vector and axial-vector con- 
tributions from dimension D > 2 operators; they contain 
an implicit suppression factor 1/M^. 

A general analysis of the relevant <5|'J'^]^ contributions 
was presented in Ref. ||^ . A detailed study of the perturba- 
tive piece S'^^^ was later performed in Ref. where a re- 
summation of higher-order corrections induced by running 
effects along the integration contour — denoted Contour- 
Improved Fixed-Order Perturbation Theory (FOPTci) — 
was achieved with renormalization-group techniques. 

The leading quark-mass corrections of dimension two 
have been studied in Refs. |,|l|,|6|,|l|,|l|,|l||; these con- 
tributions are the dominant SU(3) breaking effect, which 
generates the wanted sensitivity to the strange quark mass. 
The separate measurement of the Cabibbo-allowed and 
Cabibbo-suppressed decay widths of the r allows one to 
pin down this SU(3) breaking effect through the differ- 
ences [Ol 



T>kl 

jrjjkl _ ^T,V+A 

= 3 Sew 



\Vus 



E 

D>2 



,kl (D) 



M (D) 



(12) 



These observables vanish in the SU(3) limit, which helps to 
reduce many theoretical uncertainties. In particular they 
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s (GeV^) 

Fig. 1. Integrand of Eq. (||) with (k=0, 1=0) for the differ- 
ence ( p^ of the Cabibbo-corrected nonstrange and strange 
invariant mass spectra. The contribution from massless per- 
turbative QCD (pQCD) vanishes. To guide the eye, the solid 
line interpolates between the bins of constant 0.1 GeV'^ width. 



are free of possible (flavour-independent) instanton and/or 
renormalon contributions which could mimic dimension- 
two corrections. A detailed theoretical analysis of the lead- 
ing D = 2 and D = i contributions to ([l^) has been given 
in Ref. llEtl. 



3 Experimental data 

ALEPH has published a comprehensive study of r decay 
modes including kaons (charged, Kg and K^) p^,|30|,^, 
|3^ with up to four hadrons in the final state. Such an 
analysis is necessary in order to separate S — —1 and 
S — modes with a KK pair. 

The total branching ratio for r decay into strange final 
states is measured to be Bs = (2.87 ± 0.12)% correspond- 
ing to 

i?^,s = 0.1610 ±0.0066, (13) 

using the combined value for the electronic branching ra- 
tio, Be = (17.794 ± 0.045)%, obtained from the measure- 
ments of the leptonic branching ratios and of the r life- 
time [Q, assuming lepton universality. Since the QCD 
expectation for a massless quark is 0.1809 ± 0.0036, the 
result ( p^ is evidence for a massive strange quark. 

The strange spectral function is derived from the dis- 
tribution of the invariant hadronic mass. It is dominated 
at low mass by the K pole and the K*{892) resonance and 
at larger masses by the axial- vector resonances i4ri(1270) 
and i4ri(1400), decaying into Kmr final states. 

Fig. n| shows the weighted integrand of the lowest mo- 
ment 6E^^ from the ALEPH data, as a function of the in- 
variant mass-squared, and for which the expectation from 
perturbative QCD (pQCD) vanishes in the case of mass- 
less quarks. 



Table 1. Measured branching ratios (10"'^) for r decays into 
strange final states +Vt- The ALEPH values are from Ref. |1^ 
and the world averages (including ALEPH) froiiijlefs. 
p^ . The branching ratios for the (A'47r)~ and {K57v)~ modes 
are estimated from the measured branching ratios for the 5n 
and Gtt modes introducing Cabibbo and phase space suppres- 
sion. 



Mode 


ALEPH 


World Average 


K- 


6.96 ±0.29 


6.81 ±0.23 




4.44 ± 0.35 


4.49 ± 0.34 


K TV 


9.17 ±0.52 


8.78 ± 0.38 




0.56 ±0.25 


0.58 ±0.24 




3.27 ±0.51 


3.62 ± 0.40 




2.14 ±0.47 


2.76 ± 0.48 


K'fi 


0.29 ±0.14 


0.27 ± 0.06 


(K3n)- 


0.76 ±0.44 


0.62 ± 0.34 


7^1(1270)- -> K-LO 


0.67 ±0.21 


0.67 ±0.21 


{K4-K'y (estim.) 


0.34 ±0.34 


0.34 ± 0.34 


{K5n)~ (estim.) 


0.06 ±0.06 


0.06 ± 0.06 


Sum 


28.65 ± 1.18 


29.00 ± 1.02 



The present analysis takes into account recent branch- 
ing ratio results obtained by CLEO H) and OPAL 
thus improving the normalization of the individual contri- 
butions to the ALEPH strange spectral function. 

The largest effect is found in the Kmr final states. In 
particular, the branching ratio for the K~n~^TT~ mode, 
(2.14 ± 0.47) X 10"^ for ALEPH, becomes in the average 
(2.76 ± 0.48) X 10"'^, where the inflated error takes into 
account the poor of the fit. Also, the branching ratio 
for if^TT-Tr", (3.27 ± 0.51) x lO'^ for ALEPH, becomes 
(3.62 ± 0.40) X IQ-^ in the average. 

The resulting world average for the total strange rate 
is found to be 



Rrs^ 0.1630 ± 0.0057, 



(14) 



with a 13% improvement in precision. Table |i| summarizes 
the different contributions to Rt.s- The value for \Vus\ is 
taken from the Particle Data Group unitarity fit of the 
CKM matrix Q yielding \Vus\ = 0.2225 ± 0.0021. The 
experimental results for the spectral moments 6Rf and 
their correlations are given in Table It can be observed 
that the central values stay rather constant between A: = 
and fc = 4, while the errors decrease. This is due to the 
fact that for higher k values the spectral moments are 
relying more and more on the accurately measured K and 
i^*(892) channels. The contributions from Kmr, n > 2 
modes to the moments are negligible for k > 2. The error 
from \Vus\ is also reduced as decreases with k. The 
contributions from the various decay modes are visualized 
in Fig. 2. 



4 Phenomenological analysis 

The present analysis is performed using the general OPE 
framework including the light quark masses, the D — 
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Table 2. Measured spectral moments SR'^^ (top table: first er- 
ror is experimental, second from |Vus|) and their experimental 
correlations (bottom table). 



{k,l) 




(0,0) 


0.374 ±0.118 ±0.062 


(1,0) 


0.398 ± 0.065 ± 0.042 


(2,0) 


0.399 ± 0.044 ± 0.031 


(3,0) 


0.396 ± 0.034 ± 0.024 


(4,0) 


0.395 ± 0.028 ± 0.020 



(fc,0 


(0,0) 


(1,0) 


(2,0) 


(3,0) 


(4,0) 


(0,0) 


1 


0.94 


0.83 


0.71 


0.61 


(1,0) 




1 


0.97 


0.90 


0.82 


(2,0) 






1 


0.98 


0.94 


(3,0) 








1 


0.99 


(4,0) 










1 



4 quark mass corrections and D — 6 nonperturbative 
contributions |Q, neglecting higher dimensions. However, 
these terms are numerically insignificant compared to the 
present experimental uncertainty so that for all practical 
purposes the strange quark mass is obtained from the re- 
lation: 



(2) 

where Al., (ar) and Qkiidr) are the pQCD series, defined 
in Ref. associated with the D ^ 2 and D ^ 4 contri- 



M2 



SR. 



kl 



245" 



EW 



M4 Qkl{ar) 



(15) 



butions to SR. 



kl 



as{M^)/Tr, and 

{0\msSS~mddd\0){M^) 
-(1.5 ±0.4) X 10"^ GeV" 

(2) 



(16) 



Unfortunately, the QCD series for A^.^ turn out to 
be problematic. They exhibit bad convergence originating 
from the asymptotic behaviour of their longitudinal com- 
ponents |l^,|l^,|l9j. While the longitudinal part is known 
to third order, tne L±T series has been calculated to sec- 
ond order only. We estimate the complete third order con- 
tributions (dominated by the badly converging, but known 
longitudinal part) by assuming a geometrical growth of the 
perturbative coefficients. 



"^+^(ci+^/cf+^) 



323, of the corresponding Adler function The value 
of the strong coupling constant, as (M^) — 0.334 ± 0.022, 
is taken from the analyses |ll|,0 of the nonstrange V + A 
moments Rf. Correlations between this value of as(M^) 

(2) 

and the moments SRj^i are negligible since the uncertainty 
on the former is dominated by theory and the latter by 
the measurement of the strange component. Using these 
inputs, the pQCD series Z\^; can be displayed up to third 
order using FOPTci 0], for example: 



Pi 

CO 




s (GeVO 

Fig. 2. Integrand of Eg . (H) for the moments fc = 0, . . . , 4 and 
/ = of the difference ([I4). The widths of the interpolating 
bands corresponds to one standard deviation. 

Z\(2) ^ ;l 0390 + 0.5576 + 0.4820 + 0.4771 + . . . 



^20 (flr) = 1.1154 + 0.6432 + 0.6082 - 
^30 («r) = 1.1990 + 0.7374 + 0.7516 
Zi^2) (q^) ^ 2^ 2880 + 0.8404 + 0.9142 - 



0.6470 - 
0.8507- 
1.0928- 



(17) 



0.9734 ± 0.4811 + 0.3718 + 0.3371 



The exhibited behaviour is that of asymptotic series close 
to their point of "minimum sensitivity" and a prescription 
is needed to evaluate the expansions and to reasonably 
estimate their uncertainties. Close examination of exam- 
ples of mathematical asymptotic series suggests that a rea- 
sonable procedure is to truncate the expansion where the 
terms reach their minimum value. The precise prescription 
— cutting at the minimum, or one order before, including 
the full last term of only a fraction of it — is somewhat ar- 
bitrary and this ambiguity must be reflected by a specific 
uncertainty attached to the procedure. In this analysis we 
adopt as a rule keeping all terms up to (and including) the 

minimal one and assigning as a systematic uncertainty the 

(2) 

full value of the last term retained. It follows that the Z\^q 
series are then summed up to third order for k = 0, 1, sec- 
ond order for = 2 and first order for k = 3, 4. It can 
be remarked that the assigned truncation uncertainty is 
numerically equivalent to quoting an uncertainty of 330 
on Xg^"^, i.e., a 200% error. 

We disregard moments (fc, I) with I 7^ since they 
suffer from an increased dependence on the D > 6 non- 
perturbative terms. In addition such moments carry little 
information on m^. 

The numerical values of the D = A perturbative series 
Qki are given in Table ^. The first errors give the estimated 
theoretical uncertainties from missing higher order terms 
in the (fast converging) expansion and the second errors 
quote the changes induced by the uncertainty in the strong 
coupling. 
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Table 3. Numerical values of the relevant D — 4 perturbative 
expansions for as(M^) — 0.334 ± 0.022. The first errors give 
the estimated theoretical uncertainties from missing higher or- 
der terms in the (fast converging) series and the second errors 
quote the changes induced by the uncertainty in the strong 
coupling. 



{k,l) 


Qkiiar) 


(0,0) 


1.07 ±0.02 ±0.01 


(1,0) 


1.50 ±0.02 ±0.01 


(2,0) 


1.92 ±0.01 ±0.00 


(3,0) 


2.33 ±0.01 ±0.01 


(4,0) 


2.72 ±0.03 ±0.02 



Table 4. The strange quark mass at Mt determined from 
each of the 5R!^'^ experimental moments. The breakdown of 
the different sources of uncertainties corresponds to: experi- 
mental, \Vub\, as{M^), quark condensates, truncation of the 
perturbative series Zi^g'(aT), and renormalization scale. The 
last column gives the total theoretical uncertainty excluding 
the contribution from |Kis|. Errors have been symmetrized for 
reading convenience. 



ik,l) 


rus (MeV) 


exp. 


\Vus\ 




(T™^ (MeV) 
{msss} trunc. 


R-scale th. 


(0,0) 


132 


26 


13 


2 


4 


9 


9 14 


(1,0) 


120 


13 


9 


3 


4 


10 


11 16 


(2,0) 


117 


10 


7 


3 


6 


14 


14 21 


(3,0) 


117 


9 


8 


2 


8 


19 


16 27 


(4,0) 


103 


7 


5 


3 


9 


20 


19 29 



5 Results and discussion 

With the prescription given in the preceding section for 
the perturbative expansion Z\j,j (a^) and the data from 
ALEPH, we first derive the strange quark mass from each 
experimental SRlf moment. The results are quoted in Ta- 
ble H The values obtained for ms{M^) are rather stable 
between k = and k = A. There is however a small de- 
crease which could be of statistical nature, but could also 
indicate a deterioration of the vahdity of the OPE, since 
larger k values emphasize the low-mass contributions ren- 
dering the approach less inclusive. The breakdown of the 
error on ms{M^) into its contributions is given in Table |: 
whereas the experimental uncertainty dominates at small 
k, the theoretical uncertainty, which receives its main con- 
tribution from the truncation of the perturbative series 
(Eq. ([TtI)) and the renormalization scale, increases with 
k and dominates for fc > 1. To estimate the latter un- 
certainty, the renormalization scale is varied from 0.75Mr 
to 2Mr- All the theoretical errors are added in quadra- 
ture. Since the square of the strange quark mass is mea- 
sured (c./., H^), the errors on ms(M^) are asymmetric. 
For reading convenience, they are symmetrized through- 
out this paper, except for the final result. 

The fact that error contributions are given for the trun- 
cation and the renormalization scale, both related to the 
limited number of terms in the perturbative expansion. 



can be considered to be a conservative approach. In addi- 
tion, one can check the quoted systematic uncertainty by 
modifying the truncation procedure: cutting off the expan- 
sion one order less than the minimum term yields ms{M^) 
values ranging from 143 MeV (for fc = 0) to 127 MeV (for 
fc = 4), showing a slightly better relative stability, but de- 
viating from the nominal results given in Table by val- 
ues consistent with the quoted truncation errors. Another 
test of the handling of the perturbative series and its poor 
convergence is obtained by comparing the chosen contour- 
improved method (FOPTci) to the more standard proce- 
dure of fixed order expansion without partial resumma- 
tion (FOPT). The latter method provides ms(M^) values 
in the range 121 MeV (fc = 0) to 119 MeV (fc = 4), i.e., 
remarkably stable and consistent with the nominal values 
within the uncertainties relevant to the treatment of the 
perturbative series. 

The \Vus\ uncertainty is never dominant for any value 
of fc. Other sources of systematic effects are negligible: in 
particular, the uncertainties from Sew and higher order 
nonperturbative operators {SOq) (using the estimate of 
Ref. [0) he in the range 0.3 to 0.6 MeV. 

To optimize experimental and theoretical sensitivities, 
a combined fit is performed using several moments. The 
overall sensitivity increases up to fc = 2 while no signifi- 
cant improvement occurs for larger fc values and the fitted 
TOs(M^) value remains stable. Keeping only lower fc mo- 
ments is also justified in view of a possible breakdown of 
the OPE for higher moments and their worse perturba- 
tive convergence (c./., Eq. ([l^)). The constrained fit of 
the fc = 0, 1, 2 moments takes into account the very large 
correlations between the experimental values and yields 
the result: 

m,(Af2) = (120 ± llexp ± 8y„, ± 19th) MeV 

= (120^21) MeV, (18) 

where the theoretical error includes an additional uncer- 
tainty of (Tope = 12 MeV which accounts for the stability 
of the result when using different moments. For reading 
convenience, the errors in the first line of Eq. (^sj) have 
been symmetrized. 

The result ( p^ ) can be evolved to different scales using 
the four- loop IIGE 7- function |3^ , yielding 

rusil GeV^) = (160 1^) MeV, (19) 
TO, (4 GeV^) = (116 1^^) MeV. (20) 

As shown in the ALEPH analysis of the nonstrange r 
decays one can test the validity of the QCD analysis 
by simulating the physics of a hypothetical r lepton of 
lower mass, y/so- This is obtained by replacing by Sq 
everywhere in Eqs. (^) and (^), and correcting the latter 
for the modified kinematic factor. Under the assumption 
of quark-hadron duality, the evaluation of the observables 
as function of sq constitutes a test of the OPE approach, 
since the energy dependence of the theoretical predictions 
is determined once the parameters of the theory are fixed. 
The results of this exercise are given in Fig. ||, showing the 
variation with sq of the first four SRf moments and the 
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Fig. 3. The observables 5_R*''(so) as a function of the "r-mass" squared so confronted to the QCD predictions and the derived 
7713(80) compared to the QCD RGE running. The two bands correspond to the experimental and theoretical uncertainties. By 
construction, data and theory agree at so = M^. 



value for ms(so) derived from each moment. The bands in- 
dicate the experimental and theoretical uncertainties. The 
agreement between data and theory, perfect at sq — 
by construction, remains acceptable for lower sq values, 
down to 1.6 (2.4) GeV^ for A: = (4). For the first three 
moments used in the final determination, the running ob- 
served in data follows the RGE evolution down to about 
2 GeV^. It should be pointed out that the sq dependence of 
the theoretical prediction is obtained following the trunca- 
tion method defined in Section ^ and applied at sq = . 
If the same rule had been consistently used at each sq 
point a better agreement would have been found down to 
much lower sq values, at the price of introducing steps 
in the prediction, corresponding to dropped terms in the 
expansion (according to the prescription given in Sec. 
This observation provides another consistency test of the 
procedure. 



6 Comparison with other analyses 

Several analyses of the ALEPH strange spectral function 
have been performed in order to extract TOs(Af^). In the 



ALEPH paper |l3| , a very conservative road was followed, 
using only the L + T part since its perturbative expan- 
sion converges well. A price was paid in the experimen- 
tal sensitivity, yielding rather large uncertainties on the 
result. Also, the experimental moments were fitted, not 
only to the strange quark mass, but also to the values 
of the nonperturbative operators {50e,s). The values ob- 
tained were in reasonable agreement with our assumption 
in the present work, but unfortunately created an un- 
lucky shift to a larger value. As we have discussed 
here it is safe to neglect such contributions resulting in 
a more constraining fit. The value found by ALEPH is 
m,(M2) = (176 ^48™;; tltt ± 14„,eth) MeV, on the high 
side of the present determination for the reasons analyzed 
above. The third quoted error covers uncertainties in the 
fit and in the experimental separation of J = 0, 1 states 
necessary to work with only the L -\-T part. 

The method used in Ref. |l9| is rather close to the 
present one, with the result ms(M^) = (119 ± 12oxp ± 
lOy^^ ± 18th) MeV. Apart from using better experimen- 
tal moments, the improvements for the new analysis pre- 
sented here deals with a better treatment of the pertur- 
bative series for the different moments and an optimal 
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use of data and theory through a constrained fit of three 
moments. 

The analysis of Ref. |2| uses the ALEPH (5i?°" mo- 
ment and obtains ms{Mf) = (130 ± 27exp ± 9th) MeV. It 
advocates contour-improved resummation and employs an 
effective charge as well as effective masses absorbing the 
higher perturbative terms. The central value agrees with 
the corresponding result of the present analysis (given in 
the first line of Table The quoted theoretical uncer- 
tainty is half as small as the one derived here, but not jus- 
tified in the paper. Apparently no uncertainty is included 
from \Vus\, which should be ±13 MeV. 

Finally, the last analysis |Q using the ALEPH data 
makes use of weight functions multiplying the correlators 
m Eq. (|). These weights are tuned to improve the con- 
vergence of the perturbative series, while suppressing the 
less accurate high-mass part of the strange spectral func- 
tion. This latter feature is close to our use of higher mo- 
ments in k. Since we have observed some deterioration 
of the convergence for these moments and correspond- 
ingly a systematic shift of the derived values, it is 
not completely clear to what extent this procedure is ap- 
plicable and thus how reliable the answer is. Nevertheless 
their result is in good agreement with Eq. (p^, with a 
smaller theoretical uncertainty, apparently not including 
the effect from the arbitrary choice of the renormalization 
scale. When evolved to M^, the result of Ref. |Q reads 
m,(Af2) = (119 ± 14exp ± 12y„, ± 10th) MeV. 

Other determinations of iris have been obtained by 
analyses of the divergence of the vector and axial-vector 
current two- poi nt function correlators ||3^,p5, 36, 37|,|38|,p9|, 
||,|l|,|4|,||,g. The phenomenolo gical information on the 
associated scalar and pseudoscalar spectral functions is re- 
constructed from phase-shift resonance analyses which are 
yet incomplete over the considered mass range and need 
to be supplemented by other assumed ingredients, in par- 
ticular the description of the continuum. 

Another approach |Q considers the difference between 
isovector and hypercharge vector current correlators as re- 
lated to the 7=1 and I — Q spectral functions accessible 
in e+e^ annihilation into hadrons at low energy. A recent 
reanalysis p6[ points out the possibility of large correc- 
tions from isospin breaking, leading to significant devia- 
tions for the extracted nis value. 

Finally, lattice QCD calculations of nis are available 
(see, e.g., Ref. [0|4^ for recent reviews and references 
therein), whose results are quite spread at present. The 
average result reads Q 771^(4 GeV^) = (110 ± 25) MeV. 

The TUs determinations discussed above are compared 
in Fig. I at the scale of 1 GeV^. 

The sum of the up and down quark masses has been 
determined with Finite Energy Sum Rules |49| with the 
result (m„ -Hmd)(l GeV^) = (12.8 ±2.5) MeV. Using the 
ratio 2ms/{mu + md) = 24.4 ±1.5, obtained within 0(p'^) 
Chiral Perturbation Theory and the large Nc limit ||50| , 
this result is in nice agreement with the present determi- 
nation (|l^. 



T decays ALEPH 
[13] 
[19] 
[20] 
[22] 

this analysis 
k' SR [35,42,44] 
e*e" SR [45,46] 
Lattice [47] 
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m^(l GeV^) (MeV) 



Fig. 4. The determination of ms(l GeV^) in this work com- 
pared to the results of other analyses based on the ALEPH 
strange spectral function and on other approaches. Details are 
given in the text (SR = sum rules). 



7 Conclusions 

We have used the strange spectral function in r decays 
measured by ALEPH to extract the strange quark mass. 
The normalization has been improved by incorporating 
recent branching ratio determinations from CLEG and 
LEP experiments. The stability of the result has been 
checked by using several moments of the invariant mass 
distribution. The final value is obtained from a constrained 
fit of 3 moments, yielding 

m,(M2) = (120 ± ll,,p ± 8y„, ± 19th) MeV 
= (1201^^) MeV. 

The theoretical error accounts for uncertainties, (i), in the 
perturbative series used to 0{a^) or lower, following a 
prescription for truncating asymptotic series, (ii), in the 
OPE approach, and (in), from smaller sources. 

At the customary scales where quark masses are quoted, 
this result becomes 

m^(l GeV^) = (IGOtll) MeV, 
m^(4 GeV^) = (116^2°) MeV. 

The stability of the result is tested by varying the mass 
of a hypothetical r lepton using the measured spectral 
function. The analysis remains reliable within the given 
errors down to a mass of about 1.4 GeV. The procedure 
used for truncating the asymptotic QGD series would in 
fact improve the stability down to even smaller masses, 
with however less accuracy. 

The uncertainty on the rus result is dominated by 
theory. This does not mean that more precise data on 
Cabibbo-suppressed r decays are not necessary. In partic- 
ular the quoted error on the validity of the OPE is derived 
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from the variation of nig extracted from moments of in- 
creasing orders. It is not clear if the observed effect is 
of statistical nature and would not disappear with larger 
data samples. In general, more data, as expected from the 
B-factories presently in operation, will allow more checks 
to be performed with a possible gain in the theoretical 
uncertainty [ pT| . Thus future measurements could permit 
a more precise determination of the strange quark mass 
along the lines presented here. 
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